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Abstract 

A natural extension of the Dijkgraaf-Vafa proposal is to include 
fields in the fundamental representation of the gauge group. In this 
paper we use field theory techniques to analyze gauge theories whose 
tree level superpotential is a generic polynomial in multi-trace opera¬ 
tors constructed out of such fields. We show that the effective superpo¬ 
tential is generated by planar diagrams with at most one (generalized) 
boundary. This justifies the proposal put forward in |Q. 

We then proceed to extend the gauge theory-matrix model dual¬ 
ity to include baryonic operators. We obtain the full moduli space 
of vacua for an U{N) theory with N flavors. We also outline a pro¬ 
gram leading to a string theory justification of the gauge theory-matrix 
model correspondence with fundamental matter. 


1 Introduction 


Recently, Dijkgraaf and Vafa have proposed iBi a perturbative method 
for computing the effective glueball superpotential of several classes of = 1 
theories. The proposal instructs one to compute the planar free energy of the 
matrix model whose potential is the tree-level superpotential of the theory 
of interest. 

Being obtained via a “string theory route”, the original proposal natu¬ 
rally includes theories with helds transforming in the adjoint/bifundamental 
representation of the gauge group. One of the most natural extensions of 
this duality is to theories with helds in the fundamental representation of 
the gauge group. It turns out that this goal can be reached in a rather sim¬ 
ple way: one needs to compute only the contribution to the free energy of 
the matrix model arising from Feynman diagrams with one boundary 
Thus, the gauge theory effective superpotential is constructed as 


12 


Weff{S,A) = N,S{l 




( 1 ) 


Such a construction was successfully used to compare matrix model pre¬ 
dictions with known gauge theory results for theories with massive and mass¬ 
less havors, with Af = 1 and J\f = 2 supersymmetry |P, [T^, |T^, |27 . 

Other interesting related work has appeared in |^, |^, |^, |^, ^ 


Of a much deeper interest has been the investigation of the structure 
underlying this duality. In [Q it was proven that in a gauge theory with 
adjoint helds, the computation of the superpotential is reproduced by the 
matrix computation, diagram by diagram. Using anomaly-based arguments. 


in |T^ it was also shown that the gauge theory and the matrix model results 
are related as proposed by Dijkgraaf and Vafa. 

In this paper we prove using held theory techniques that the ehective 
superpotential of a theory with a tree level superpotential being a generic 
polynomial in multi-trace and baryonic operators is generated by planar di¬ 
agrams with at most one (generalized) boundary. This justihes the proposal 
put forward in |]T^. We then proceed to extend the gauge theory-matrix 


model duality to include baryonic operators Q We obtain the full moduli 

^As we were preparing the manuscript, appeared which has some overlap with our 
discussion of baryons. 
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space of vacua for an SU{N) theory with N flavors. We then outline a pro¬ 
gram for the string theory justihcation of the extensions of the Dijkgraaf-Vafa 
duality. 


2 Field theory analysis of the effective super¬ 
potential 

In this section we analyze in detail the effective superpotential as a function 
of the glueball superheld S = Tr{W‘^) and the various coupling constants 
that exist in the theory, and show that in theories with helds transforming 
in the fundamental representation of the gauge group the superpotential 
is generated entirely by Feynman diagrams with a single boundary. The 
analysis is similar to the one described in [^. We will hrst consider tree 
level superpotentials built out of traces of products of quark bilinears. This 
analysis proves the proposal put forward in |^, that only diagrams with 
one boundary contribute to the effective superpotential. We will then add 
multi-trace deformations as well as baryonic operators. 

We begin by considering a theory with a fairly arbitrary, polynomial, tree 
level superpotential: 

= , ( 2 ) 

1>1 

where the trace is over the flavor indices. If no tree level superpotential is 
present, the global symmetry of the theory is SU {Nf) x SU {Nf) x U{l)iX 

U{l)iXU (1)r. The introduction of the superpotential above, with scalar 
couplings, breaks this symmetry group to 


SU{Nf) X U{1), X U{1), X U{1)r 


( 3 ) 


As in the case of a theory with adjoint helds [0, these symmetries are 
not sufficient to completely determine the effective superpotential. This is 
easy to see by considering the following rather unusual charge assignment for 
the various helds and couplings: 



SU{Nf) 

[/(I) 

[/(I) 

mn 

A 

g: 

Nf 

1 

0 

1 

1 

Q: 

1 

0 

1 

1 

1 

9i ■ 

1 

-1 

-1 

2-21 

3-21 

S : 

1 

0 

0 

2 

3 


( 4 ) 
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The last column in the table above represents the engineering dimension of 
the corresponding field. 

A simple counting argument shows that all higher loopQ Feynman dia¬ 
grams with insertions of the glueball superheld are hnite. We are thus look¬ 
ing for a function which depends only on S and gi. Furthermore, since this 
function is generated perturbatively, neither one of its arguments is allowed 
to appear at a fractional power. It is then not hard to see that the basic 
invariant combinations are 


(V) l>2 (5) 

Qi 

and they also have vanishing scaling dimension. Thus, the most general 
effective superpotential which can be generated is 

«'ee = SF{^) ( 6 ) 

Ul 

where F is an arbitrary analytic function of its argument^. A generic term 
in the series expansion of Weff is 




n 

i=l 9l 


(7) 


where the li i = 1,... ,n are not necessarily distinct and are numerical 
coefficients. It is clear from the coupling constant dependence that this term 
is generated by a Feynman diagram with n vertices with the corresponding 
couplings gl^ i = l,...,n. The F-dependence of this term also hxes the 
number of index loops in the diagram. Indeed, since S involves a trace, it 
is clear that each index loop will generate at most one factor of the glueball 
superheld. This implies that the diagram will have at least 

n 

i = l + Eft-l) (8) 

i=l 

^The one loop graphs are the only divergent ones; their contribution can be taken into 
account using anomaly arguments ||l5|| . An attempt in this direction was made in 

^The appearance of inverse powers of gi is allowed perturbatively, as it can be accom¬ 
modate by rescaling the quark fields. Indeed, absorbing in both Q and Q implies that 
the relevant couplings to consider are ^. 
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index loops. As in the case of diagrams with helds transforming in the adjoint 
representation of the gauge group, it is possible to show that, for a Feynman 
diagram with vertices generated by (|^), which has genus g and b boundaries, 
the number of gauge index loops is given by 

n 

Lg = 2-2g-b + Y.ik-l) . (9) 

i=l 


Thus, since the theory we are interested in has no helds in the adjoint rep¬ 
resentation (and thus each diagram has at least one boundary), the only 
way for the relevant Feynman diagrams to have the number of index loops 
required by is to be planar and have exactly one boundary. This proves 
the proposal put forward in . 


It is easy to add multi-trace deformations to the tree level superpotential 
in equation (Q). The coupling constants now carry several indices, counting 
the number of traces as well as the number of quark bilinears in each trace. 
Thus, the terms with I quark bilinears are: 


y 9S?„,„rr[(Q0)"].--rr|(QQ)””| (10) 

Sl-\- 


It is clear that the charges of gs}...sm under the symmetry group are the 
same as the charges of gi in the table above. Thus, the generalization of 
to this case is 

q{i) 

- (V) / > 2 , (V) m . (11) 

9i 

As in the previous case, the Feynman diagrams generating the term propor¬ 
tional to must have 1 -|- J^iih — 1) index gauge loops. Unfortu¬ 

nately, the relation between the type of vertices used in the diagram, the 
number of gauge index loops and the number of boundaries becomes more 
complicated than (^). It is however not hard to see that the number of 
index loops remains unchanged if one replaces a single-trace vertex with a 
multiple-trace vertex while keeping the genus unchanged. 

Graphically, this corresponds to a degeneration of a disk by pinching its 
boundary, as is shown in hgure |^. On the left-hand-side of that hgure we show 
a vacuum diagram constructed out of a 4-point and a 6-point single-trace 
vertices. The outer line represents the boundary (corresponding to a trace 
over flavor indices) while the inner lines represent gauge index loops. On the 
right-hand-side we show two possible degenerations of this disk. Following 
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route (a), we replace the single-trace 4-point vertex by a double-trace one. 
Following route {b), we replace the single-trace 6-point vertex by a double¬ 
trace one. Since these surfaces lie at the boundary of the moduli space of 



Figure 1: Degenerations of the disk. 


the disk with holes, it is clear from this perspective that the latter diagrams 
should contribute on the same footing as the former ones. With this dehnition 
of boundary components, it follows that equation (||) still applies and the 
diagrams contributing to the superpotential have exactly one boundary. 

Let us now turn to including baryonic operators 


if Nf = Nc and 


B = detQ B = detQ (12) 


B^ = 

O 



Bi = 

■■^Ncnh 

V ai • 


( 13 ) 


if = W + 1. We will repeat the analysis in the beginning of this section 
and construct the basic invariants which can be built out of these operators. 
We will deform the initial theory only by terms linear in the baryons. More 
general deformations can be treated similarly. The charges of the baryon 
sources are then: 


b: 

SU{Nf) 

1 

m 

-iVc 

m 

0 

mR 
2-W 

A 

3-A, 

(14) 

b: 

1 

0 

-iVc 

2-W 

3-A, 



It is easy to see that the basic invariants are 

b b 




gi 


N ^/2 N ^/2 


9i Igr 9i 


^{l-l)+n{Nc- 2 ) 


( 16 ) 
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where n is an arbitrary number. Thus, the superpotential is a series expan¬ 
sion in Ii^n- We can then immediately see the Feynman diagram origin of a 
generic term in this series: 

^ ^ cjVc-2 

N,/2 N,/2 ^ 

i/1 i/l 




If all rii vanish, this term matches equation (|^) and each color index loop 
contributes one power of S. Consider then adding a baryon and anti-baryon 
to the diagram giving this contribution. To produce a IPI diagram we need 
to break one of the initial color index loops and connect it to one Q and one 
Q from the baryons. Thus, we hnd that the number of index loops increases 
with Nc — 2 ii the quark lines are connected in a planar way. Adding more 
baryon-anti-baryon pairs follows the same pattern and recovers the second 
factor in equation ([l6|) . 

We would like to point out that the notion of planarity we use here 
is slightly more general than usual. Due to the presence of the Levi-Civita 
symbols at the baryon vertices one is allowed to use their symmetry properties 
to reorganize the internal lines of a diagram. 

With hindsight, this result is not surprising in light of the discussion 
concerning multi-trace deformations of the initial theory. Since baryon-anti- 
baryon pairs can be re-expressed as sums of multi-trace operators, it follows 
that the counting of powers of S for correlation functions involving insertions 
of such pairs is identical to the counting of powers of S for correlation func¬ 
tions involving insertions of multi-trace operators. There we saw that each 
color index loop gave a single S and the diagram was required to be planar 
and have a single boundary in the generalized sense described above. 

We have thus proven that in theories with a tree-level superpotential 
containing multi-trace operators and/or baryon operators, the effective su¬ 
perpotential as a function of the glueball superheld is generated by planar 
Feynman diagrams with a single boundary in the generalized sense dehned 
above. 


3 On the dynamics of massless fields 

One of the most puzzling aspects of the Dijkgraaf-Vafa proposal relating a 
gauge theory with a matrix model is the fact that nonperturbative physics in 
the former can be obtained by doing perturbative computations in the latter. 
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To relate the two theories, we are instructed to hrst compute the free 
energy of the matrix model and construct from it the gauge theory effective 
superpotential. As it is the case with all dualities, there exists a dictionary 
between objects on both sides, but these objects are not the same. One 
can only compare correlation functions involving these objects, but not the 
objects themselves. Therefore, to make contact between the two theories, one 
needs to integrate out all helds and to compare the effective superpotentials 
as (complicated) functions of the coefficients of the tree level superpotential 
and the dynamically generated scale. 

It seems therefore that the only way one can compute gauge theory super¬ 
potentials is to use the matrix model to compute the effective superpotential 
in the “deepest infrared” in the Wilsonian sense, and then to integrate in 
all these helds. This seems to limit the use of matrix models in describing 
theories which one cannot obtain by integrating in helds. 

This procedure of reconstructing the gauge theory superpotential out of 
the matrix model results may seem somewhat cumbersome. Recently, for 
theories with tree-level superpotential expressed only in terms of quark bilin- 
ears (mesons), a method was proposed [|T^ to link the matrix model and the 
gauge theory directly, bypassing the detour through the “deepest infrared”. 
The method consists of introducing in the matrix path integral a S function 
which relates the matrix model quark bilinear with the gauge theory meson. 
While this method seems rather ad hoc, as it identihes two quantities which 
have nothing to do with each other (they are related by a dictionary, but 
they cannot, in principle, be identihed), it implies that the full effective su¬ 
perpotential is given by the sum between the tree level superpotential and an 
Affieck-Dine-Seiberg nonperturbatively-generated term. The fact that this is 
the correct answer can be argued using standard gauge theory arguments. 

It is unclear at this time how to extend this program to directly include 
superpotentials with terms depending of several types of helds (e.g. Yukawa 
coupling), with terms which cannot be written in terms of quark bilinears 
(e.g. baryon operators in SQCD with Nf > Nc). 

This procedure is moreover not unique. In the appendix we present an 
alternative procedure relating the matrix model and the gauge theory which 
also gives the correct gauge theory result and moreover treats different rep¬ 
resentations on the same footing. Both the procedure of |T^ as well as the 


one described in the appendix are best thought of as shortcuts which relate 
the two sides before integrating out the helds on either side. It would be 
quite interesting if one can justify them to be more than this. 


7 






Let us proceed by reviewing the method proposal of |^. The idea is to 
identify the matrix model gauge invariant quark bilinear with the held theory 
meson, by introducing a constraint for the quarks in the matrix integral. 
More formally, 

= J DQ,DQ, 6 iM,, - ^ ( 17 ) 


where the 6 function enforces a direct relation between the high energy de¬ 
scription (in terms of the independent helds Qj" and the one in 

terms of composite helds Mij. 


Because of this 6 function, any tree level superpotential which can be 
expressed as a function of the mesons automatically gets out of the integral. 
The remainder of the integral gives the volume of the constrained quark in¬ 
tegration domain. As in the original DV prescription, the glueball superheld 
is identihed with the dimension of the matrix model matrices (at least in 
the case of U (M) matrix model gauge group) after it was taken to be large. 
Adding the usual Veneziano-Yankielowicz terms (to account for the gauge 
held dynamics |T^) and integrating out the glueball superheld one obtains 
the AHleck-Dine-Seiberg superpotential. The net result of (^) is to always 
give the meson ehective potential 


^eff — Wads + Wtree ■ ( 18 ) 

For simple tree level superpotentials this is indeed what one expects to 
hnd in gauge theory [^. However, it is not a priori clear that for a generic 
tree level superpotential the ADS superpotential is the only possible nonper- 
turbative contribution. If this were not the case, the procedure (^) would 

fail. 

We can outline the sketch of a gauge theoretic argument that indeed there 
is no nonperturbative contribution besides the ADS piece. Let us assume a 
Wtree of the form 

Wtree = Y.9^,kTr[Xf ( 19 ) 

i 

Besides the ADS superpotential, one can expect functions of the dimension¬ 
less and chargeless ratios of coupling constants (e.g. 32311 etc.) to 

appear multiplying Wads- However, these functions need to be expandable 
in positive integer powers of the coupling constants around the origin. More¬ 
over they need to be holomorphic in all couplings. Hence, these functions 





would be holomorphic functions which go to zero at oo, and would thus be 
zero everywhere. Thus, the ADS superpotential will be the only nonpertur- 
bative contribution to the superpotentiah 

As pointed out in the beginning, the procedure (^) is an effective way to 
bypass integrating out on both sides when the potential can be expressed in 
terms of mesons only. However, extending it to helds in other representations 
appears challenging. 

For example, in theories containing massless adjoints the number of mesons 
one can construct becomes very large |^, |^. Multiple mesons can appear 
also in the case one considers theories with massless quarks in the symmetric 
or the antisymmetric representations of the gauge group. In order to capture 
the dynamics of these mesons, one would need to add them by hand in (|^. 
This would further constrain the quark integration volume. As before, any 
superpotential depending on the mesons would get out of the integral, and 
the constrained volume would give the nonperturbative superpotential as a 
function of all these mesons. 

Nevertheless, if one tries to apply this procedure to baryons there is no 
way to get sensible results. Indeed, a gauge theory baryon is a combination 
of Nc quarks. However, the matrix model prescription instructs us to take 
the matrix model gauge group to have large rank, unrelated to the gauge 
theory Nc- Thus, depending upon how one chooses to treat it, the naive 
baryon operator either vanishes or is not gauge invariant. A matrix theory 
baryon should be composed of M quarks, where M is the dimension of the 
matrix. It is quite clear that the matrix model baryon and the gauge theory 
one have nothing in common with each other. 

One can actually see that baryons are a problem not only for various 
computational “shortcuts”, but also for formulating the correspondence in 
the picture in which one hrst goes to the deep infrared of both theories. One 
should add a term of the form in the gauge theory 

and integrate it out. Nevertheless, it is quite clear that such a term does not 
make any sense in the large M U{M) matrix model where the quarks trans¬ 
form in the fundamental representation of U (M) but their number remains 
hxed to Nf. The next section is devoted to solving this problem. 
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4 Baryons 


In this section we try to understand how baryons enter in the gauge theory 
- matrix model correspondence. As we explained above, introducing baryon 
sources both in the large M matrix model potential and in the gauge the¬ 
ory superpotential cannot be done. Thus, one has to take a step back and 
consider the motivation behind the DV proposal. 

The proposal was hrst formulated via a string theory route. Unfortu¬ 
nately, understanding the role of fundamentals from this standpoint is rather 
challenging. In the last section of this paper we will return to the question of 
how one can understand theories with helds in the adjoint and fundamental 
representations via geometrical transition. 

The purely gauge theoretical motivation of this duality was formulated 
in [^], and in a nutshell can be summarized as: “planar diagrams in theories 
with adjoints compute holomorphic quantities.” 

If one adopts this outlook, the matrix model is merely a tool to system¬ 
atically compute planar diagrams. Phrased in this way, the “dictionary” 
between the two sides does nothing but select the planar diagrams in the 
gauge theory. Of course, it is still amazing that planar diagrams compute 
quantities which one normally obtains nonperturbatively, but the “miracle” 
now pertains only to the gauge theory side. 

The extension to theories with fundamentals and no baryons was conjec¬ 
tured in to be: “planar diagrams and diagrams with one quark boundary 
compute holomorphic quantities”. In section 2 we gave a gauge theoretical 
proof, along the lines of |T^, that indeed only diagrams with one boundary 


contribute to the low energy effective superpotential. As a side remark, we 
should note that there is still no string theory motivation for this extension. 

In the same spirit as 0, the matrix model with quarks and the correct 
dictionary between the two sides |^, should be thought of only as a tool 
to systematize the gauge theory computation of these diagrams with bound¬ 
aries. 

If one introduces now diagrams with baryons, it is clear that they are 
quite different from diagrams with quark loops. Indeed, the quark lines are 
not closed, but rather go from the baryons to the antibaryons. Nevertheless, 
we have shown in section 2 that a pair of baryons can be thought of as 
a multi-trace operator. We have also introduced a generalized notion of 
boundary for this kind of operators, and we have proven that only planar 
Feynman diagrams with one generalized boundary can generate the effective 
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superpotential. 

The aim is therefore to analyze in a systematic way these planar dia¬ 
grams. Unfortunately, a large M matrix model appears not to be any more 
the right “tool” for selecting these diagrams. Baryon terms of the form 


T-Nc 


cannot be introduced in the Matrix Lagrangian, 
because the matrix quarks are 1 x M matrices, and the combinations above 
are not matrix scalars. 

Thus, the only types of gauge theory superpotentials one can insert in 
the large M Matrix Lagrangian are the ones which can be expressed in terms 
of color singlet quark bilinears. As we have seen in the previous section, in 
the this case the shortcut of |]T^ can be applied, and one simply obtains the 
sum of the ADS superpotential and the tree level one. 

Since the large M matrix model cannot be used as a tool to select planar 
Feynman diagrams when baryons are present, the only way to count these 
diagrams is to select them by hand in a matrix model whose gauge group 
U{M) has the same rank as the gauge theory gauge group U{Nc). We stress 
that this complication is necessary only when baryons (or more generally 
operators which exist only for fixed values of the N^) are present. 

By identifying again g'^M with S, one obtains an effective superpotential 
Weff(S,mi,b,b...). Counting these diagrams is quite nontrivial, especially 
for baryons with many quarks, or for theories containing other terms besides 
quark masses and baryons. 

If one integrates S out of the matrix result Weff{S,mi,b,b...), and inte¬ 
grates in the quarks one should obtain the full gauge theory result. 

It is also possible to see that, due to charge conservation, a vacuum 
diagram contains an equal number of baryons and anti-baryons. The product 
of a baryon and an antibaryon can be reexpressed as a multiple trace operator, 
by using the fact that a product of e tensors can be expressed as a sum of 
products of Kronecker S functions. Therefore, any nonzero vacuum Feynman 
diagram containing baryons can be written as a diagram containing multi¬ 
trace operators, and as such the large M limit gives sensible results. 

Alternatively, one can integrate out a quark linking the two baryon, and 
obtain an effective di-baryon operator which again makes sense in the large 
M limit. 

It is not hard to see that this discussion combined with the fact that multi¬ 
trace tree level potentials in the gauge theory only receive Wads as a non- 
perturbative correction extend the correspondence to theories with baryons. 

In the next section we will compute the effective superpotential for an 


11 



U{2) gauge theory with 2 flavors, and show that one recovers the expected 
gauge theory result. We then use the shortcut (^) and recover the quantum 
moduli space of SU{N) SQCD with N flavors. 


5 Effective superpotentials for mesons and 
baryons from Matrix Models 

After discussing the mesons and the baryons, we now turn to calculating 
the gauge theory superpotential from matrix model. We start with a brief 
review of the held theory results for theories where M = 2 SUSY is broken 
to W = 1. 


5.1 Field Theory results 

Consider^A/” = 1 SU{Nc) theories with Nf quarks in the fundamental repre¬ 
sentation, obtained from M = 2 SU{Nc) by turning on a mass term for the 
adjoint chiral multiplet: 

W = V2Tt(Q(^Q) + ^Tr02 (20) 


To discuss the moduli space of this theory, we distinguish between the case 
when the mass of the adjoint held is small (the adjoint held remains as a 
held in the A/” = 1 theory) and the case when the mass of the adjoint held is 
large and the adjoint is integrated out and we obtain M = 1 supersymmetric 
QCD with the dynamical scale: 


A 


3N^-Nf 

N=l 



2N^-Nf 

N=2 


( 21 ) 


After integrating out the adjoint held, one obtains: 


W 


Tr[(W)(QQ)] 


iTr(OQ)Tr((3(3) 


( 22 ) 


For cx), this superpotential approaches zero and we then obtain Af = 1 
SQCD. 

"^Unlike the rest of the paper, in this section we consider the gauge group to be SU (W) 
rather than U{Nc). We do this to conform with the papers referred to in this section. All 
relevant results also hold for U {Nc) ■ 
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In a detailed discussion was presented, relating the vacua of 

A/” = 2 SUSY and M = 1 SUSY theories. Without going into details, we 
recall now the main results regarding the M = 1 vacua: 

• for Nf < W — 1, the strong gauge dynamics corrects the tree level 
superpotential by the addition of the ADS superpotential @). The moduli 
space is the extrema of this superpotential and, at hnite values for the mass 
of the adjoint held, it has been shown to have 2 types of solutions 
when all the diagonal entries of M = QQ are equal to 


35 


one 


Mu = {- 


•Nr-N., 


''ffiA 


(23) 


/ 


and the other one when there are two different diagonal entries (there are 
only two possible values for the diagonal entries). The presence of this two 
choices remains valid if we add masses for the quarks and integrate them out. 
The hrst case has been discussed in matrix model in |^, . The second one 

can be treated similarly. 

In a pure N = 1 SQCD, with vanishing tree level superpotential, all these 
vacua run away to inhnity. 

• for Nf = Nc, the situations changes because the classical moduli space 
of vacua is not only parametrized by mesons, but also by baryons. In the 
case Nf = Nc, the classical moduli space of vacua is modified quantum 
mechanically by 


detM - BB = (24) 

In this case there is no nonperturbatively-generated correction to the effective 
superpotential; the constraint (^H) is implemented with a Lagrange multiplier 
X. The exact effective superpotential is then 


Weff = W = X(detM -BB- A^^^) + 


Tr(M^) - —Tr(M)Tr(M) 
L Nq 


.(25) 


As before, M can again have at most two different eigenvalues; for A = 0 
we return to the situation discussed before, for X ^ 0 one has new branches 
which did not exist for Nf < — 1: 


Mu = A 


2Nc 


( 26 ) 


for BB = 0 or 


MY -BB = a 




(27) 
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for BB 7^ 0 

• for + 1, the classical moduli space is not modified, but the 

number of baryons increases. In this range, there is a dual magnetic descrip¬ 
tion based on the gauge group SU{Nf — Nc) with Nf fundamental flavors and 
gauge invariant helds M, the Seiberg-dual theory Discussions on the 
realization of Seiberg’s duality for mass deformed theories in the context of 
gauge theory-matrix models appeared in 
a duality for massless theories. 


17 , but it is unclear how to obtain 


5.2 Quantum Moduli Spaces from Matrix Models 

As explained earlier, the helds appearing in the matrix model are not directly 
related to the dynamical gauge theory helds. Thus, one cannot directly 
identify the matrix model ehective action in which only some helds were 
integrated out with the gauge theory superpotential. Instead, to hnd the 
gauge theory ehective superpotential using the matrix model, we take the 
following steps: 

1) First we deform the gauge theory by including mass terms for all helds 
as well as sources for the operators describing the directions in the moduli 
space we are interested in. In the case of mesons these two coincide, as the 
mass of a quark acts as a source for a meson as well. 

2) We compute the planar free energy of the matrix model associated to 
this deformed theory. 

3) We equate the resulting function with the gauge theory ehective su¬ 
perpotential at energy scales below the lightest mass present in the deformed 
theory. The ’t Hooft coupling of the matrix model is identihed with the 
gauge theory glueball superheld. Furthermore, the Veneziano-Yankielowicz 
superpotential is added to take into account the gauge dynamics. 

4) We integrate in gauge theory helds by appropriately performing Leg¬ 
endre transforms 

5) We restore the masses of the original helds, by appropriately taking 
the deformation parameters to zero. 

To illustrate this program let us discuss in detail an SU{2) gauge theory 
with 2 havors, Qi and Qi with i = 1, 2. This means that the adjoint held 
(appearing in the discussion in the previous subsection) has been integrated 
out and we work within M = 1 SQCD. We return at the end of this section 
to the case when the adjoint held has hnite mass. We start with a massless 
theory whose tree level superpotential is given by the baryon and antibaryon 
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sources: 


tree 


d^e 


bQi°‘Q2^€al3 + bQi^Q2^^ap 


+ he 


(28) 


The first step then instructs us to dehne this theory as the massless limit of 
the theory with the superpotential above deformed by mass terms: 


Weleff = / d^9 


+ m(52“Qa] + he 


(29) 


As discussed earlier, it is clear that Wtree cannot be used as part of the 
potential of a large N matrix model (or any other theory for that mater), 
because it manifestly vanishes. Thus, the way to proceed is to consider an 
SU (2) matrix model and compute only the planar free energy. Thus, we have 
the following matrix model potential: 


WM = {Ql",Q2^,Qa\Qa) 


( 0 heap rn 

heap 0 0 

m 0 0 

Vo m be^h 


m 


\ 

/g/\ 


Q2h 


Qp^ 

/ 

w// 


= CpMQ 


(30) 


Since this potential is quadratic in the dynamical helds Q and Q, the free 
energy is trivial to compute. The result is: 




-InZ 


In 


■ 1 
^2x2 


1/2 


detM 


2 X 2 X - 
2 


ln( 


vn? — bb 
A2 


(31) 


In the expression above the hrst factor of 2 is due to having two flavors while 
the second factor of 2 is due to having two colors. This latter numerical 
coefficient combines with the coupling constant into the glueball superheld. 
Therefore the gauge theory superpotential is: 


Wg = N,S 


1 — In 


S 

A3 


1 AT ^-bb + 

+ ' 


(32) 


where in the case at hand Nf = Nc = 2 but we wrote them explicitly to 
stress the functional dependence, and we inserted A for dimensional reasons. 

Now we perform the Legendre transform required for integrating in the 
quarks in the form of the gauge theory mesons Xi = QiQ^ and baryons 
B = Qi^Qh^ap and B = Qi^Q^^ap- This implies that we deform Wg by 


Wl = -mNfX -bB-bB 


(33) 
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and integrate ont the sonrces m, b and b. Here we also considered all mesons 
to be identical. Eliminating the sonrces by their eqnations of motion and 
nsing the eqnality between Nc and Nf we hnd: 

Weff(X,B,B,S) = 2S\n -=- . (34) 

^ ^ J -BB + X^ 

Upon eliminating S by its eqnation of motion we hnd the correct gange 
theory snperpotential 

W,ff = 0 (35) 

together with the qnantnm corrected modnli space: 

det X -BB = A^ (36) 

where we wrote in the SU{2) x SU{2) invariant form. This reprodnces 
the gange theory resnlt. 

5.3 The Effective Potential from Integrating Out the 
Massive Quarks 

Let ns now consider the general case of an SU{Nc) theory with Nf = Nc 
havors and discnss the relation between the matrix model and gange theory. 
Since we are interested in incorporating baryonic operators in the matrix 
model, we deform the gange theory action by sonrce terms for them. Then, 
the idea we will pnrsne is the following. Using matrix model reasoning we 
rednce the matrix model with Nf = Nc flavors to one with Nf = Nc — 1 
flavors. Then, for this new matrix model we constrnct the corresponding 
gange theory snperpotential. Then we integrate in, as gange theory qnark, 
the qnark which was integrated ont as a matrix model one as well as the 
baryons. To make a long story short, we will recover the known gange theory 
resnlt, i.e. a vanishing effective snperpotential and a qnantnm-deformed 
modnli space. 

To begin with, the deformed tree level snperpotential is: 

Wdef = + bB + bB (37) 

where the baryon operators B and B are dehned in eqnation (|T^). Since B 
and B are antisymmetric prodncts of the Nf qnarks, we can split them into 
prodncts of the massive qnark and the remaining massless qnarks: 

B = (38) 
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where B°‘ is product oi Q 2 i -'' tQn^- The same applies for the antibaryon: 


B = (39) 

With these slight adjustments, the superpotential reads 

Wdef = mQiQi + (40) 

Since the massive quarks appear only linearly and quadratically, it is easy to 
compute the matrix model effective action by integrating out Qi and Qi as 
matrix model helds. The computation is entirely straightforward and yields 

771 ho ^ ~ 

where we identihed the size of matrices in the matrix model with the gauge 
theory glueball superheld. 

Because i?" and S" are products of W — 1 flavors, 


B^Ba — det(M7v-=Vc-i) 


( 42 ) 


To obtain the gauge theory superpotential one should also integrate over 
the massless quarks. Instead of computing the free energy with this effective 
potential and then integrating in the gauge theory quarks, we can just apply 


the prescription 
gives: 


which we reviewed in section 0. This straightforwardly 


Weff{S,MM^-um,b,b) 


S (-In^ + 1) - ^ In 

^ A3 ^ 




det{MNf=Nc-i) 

\2Nc-l 


(43) 


The term proportional of the logarithm of the mass of the quark which was 
integrated out combines with the initial dynamical scale A to give the scale 
of the theory with Nf = Nc — 1 quarks, just as in a purely held-theoretic 
analysis. Then, to construct the gauge theory superpotential in terms of 
meson helds only, we integrate out the glueball superheld and hnd the usual 
ADS superpotential: 




A 2 Nc +1 

^^Nf=Nc-l 


det{MN^=Nc-i) 


— {B^Bc 
Am ^ 


(44) 
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The last step is to integrate in the gauge theory quark corresponding to 
the matrix model quark which was integrated out, and restore the last meson 
as well as the two baryons. The steps are identical to those in and we 
will not repeat them here. The conclusion of this short computation is that 
we recover a vanishing superpotential as well as the usual modification of the 
classical moduli space of the theory. 


6 Geometric Transitions, Matrix Models and 
Gauge Theories for Fundamental Quarks 

The geometric transition was first formulated in type IIA string theory from 
one side with D6 branes wrapped on the cycle of a deformed conifold 
to another side with fluxes on the cycle of a resolved conifold (see 
l3^ , |39| , ^ for subsequent work on type IIA). 

The mirror picture is a type IIB transition from D5 branes wrapped on 
the of a resolved conifold to fluxes on a cycle of the deformed conifold. 
The type IIB picture was easier to generalize, by using Af = 1 deformations 
of resolved M = 2 ADE singularities |^, ^ . 

The geometrical picture was the backbone of the recent exciting devel¬ 
opments in matrix model/gauge theories. In [|l|, the model of has been 
used to relate matrix model variables with geometrical variables and with 
gauge theory variables. In the comparison between geometry, matrix models 
and gauge theories, what needs to be compared is: 

• in field theory we consider the quantum moduli space. One of the tools 
to study this are the Seiberg-Witten curves. Going beyond the field theory, 
gravitational corrections to field theory have been discussed in the matrix 


model description [47, 4S 


• in geometry, we consider the 1-complex dimensional curve obtained by 
reducing the deformed geometry on an 3“^ fiber. 

• in the matrix model, we consider the equation for the resolvent of the 
matrix. 

If we are to compare just field theory and matrix model, in Q has 
been suggested that the planar diagrams in matrix models give the non- 
perturbative corrections in the field theory (we have discussed extensively in 
our paper about this). In the line of [^, it is interesting to get a connection 
between field theory and matrix model, through geometrical quantities. It 
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would be interesting to also describe matter (besides bifundamental fields) 
in the geometrical setting. It is not yet known how to go through the tran¬ 
sition with fundamental matter, and in the present section we will outline a 
possibility, leaving the details for a future work. 

To begin with, we have to know how to include fundamental matter in 
the geometrical set-up. The fundamental matter would correspond to ei¬ 
ther additional D5 branes wrapped on additional non-compact cycles of the 
conifold or to additional D7 branes which are partially wrapped on the 
cycles. In order to see this, one can use the MQCD picture for the geomet¬ 
rical transition |50| , Q. In the absence of the fundamental quarks, 
the conhguration with D5 branes wrapped on cycles goes into a brane 
conhguration with D4 branes between two orthogonal NS branes [Q (there 


are two orthogonal lines of singularity in the geometry, which become two 
orthogonal NS branes after the T-duality). In the brane conhguration, the 
fundamental matter can be introduced as either semi-inhnite D4 branes end¬ 
ing on one of the NS branes or as D 6 branes in between the NS branes. The 
D4 branes are better suited to describe massive havors and the D 6 branes 
are better suited to describe massless havors. In the geometric picture the 
semi-inhnite D4 branes become D5 branes on a non-compact P^ which is 
distanced from the exceptional P^ by a distance equal to the mass of the 
havor [Q. We do not yet have a clear picture on the geometry side for the 
D 6 brane, but we can use the geometry/brane conhguration duality in order 
to discuss the transition. 

For massive matter, in the infrared we go to the ehective theory where all 
the massive matter has been integrated out and the sizes of the diherent 
cycles on the geometric side are related to the diherent gluino condensates 
in held theory, which are in turn proportional to the masses of the quarks. 
If one of the masses of the quarks becomes zero, the size of the cycle 
becomes zero and the geometric transition fails. In this case, we should use 
the brane conhguration with D 6 branes instead of the one with semi-inhnite 
D4 branes. By lifting it to MQCD and going through the transition, the 
moduli space of the held theory should be recovered from the geometry. As 
the transition consists on closing the interval between the two NS branes 
(going from a curved M5 brane to a plane M5 brane |^), we expect to 
get a geometry described by an equation xy = a m. the (x, y) plane. But 
section 5.2 of tells us that the only case when xy 7 ^ 0 is when the mass 
of the M = 2 adjoint held is not inhnite. This is consistent with the fact 
that there is no supersymmetric vacua for the corresponding M = 1 theory. 
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In the case of infinite mass for the adjoint held, after integrating it out in 
the matrix theory, we remain with decoupled integrals over the massless 
fundamental quarks (with no tree level superpotential), and this gives the 
ADS superpotential which removes any supersymmetric vacua. 

In order to obtain a vacuum, we have to keep the mass of the M = 2 
adjoint held as a hnite quantity. In this case the matrix integral becomes 
harder to compute, as there are extra couplings between helds. This is cur¬ 
rently under investigation. 

We can ask whether there exists a geometrical picture corresponding to 
hnite mass for the M = 2 adjoint held. Remember that the starting point of 
the geometrical discussion was the case when the mass of the M = 2 adjoint 
held is cx). In the brane conhguration this comes from the fact that the NS 
branes are orthogonal and this comes from the orthogonality of two lines 
of singularity in the geometry. As we know from brane conhgurations, the 
angle between the NS branes is related to the mass of the adjoint, therefore 
a hnite mass for the adjoint held can be obtained only when the NS branes 
are not orthogonal which would imply a geometry with non-orthogonal lines 
of singularity. 

7 Conclusions 

In this paper we discussed several aspects of the Dijkgraaf-Vafa relation be¬ 
tween matrix models and gauge theories. Our investigations focused on the 
helds transforming in the fundamental representation of the gauge group. 
Using held theory arguments (based on symmetries and holomorphy) we 
proved that the ehective superpotential for the glueball superheld is gen¬ 
erated on the gauge theory side by planar (in a slightly generalized sense) 
Feynman diagrams with exactly one boundary. Our result holds for a tree 
level superpotential given by a generic polynomial in multi-trace operators 
and baryonic couplings. It is then likely that, using techniques similar to 
those in |^, the gauge theory computation be explicitly reduced to a matrix 
model computation. A slight subtlety in the diagrammatic analysis is due to 
the existence of graphs which are not obviously planar, but can be mapped 
into manifestly planar ones using the symmetries of the baryonic vertices. 

We also analyzed the inclusion of baryonic sources in the matrix model 
computation. This solution turned out to be rather subtle, as one cannot 
include the baryon operators in a large N matrix model computation. In- 
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stead one must keep the rank of the matrix model gauge group equal to the 
rank of the gauge theory one. Then, the restriction to planar diagrams is not 
automatic any more and must be done by hand. We explicitly performed the 
computations for U{2) SQCD with two flavors and then used a “computa¬ 
tional shortcut” to analyze the general case of U{N) SQCD with N flavors. 
We recovered the known gauge results for the quantum moduli space. 

As the geometry was a crucial tool in deriving the relation between matrix 
models and gauge theories in |l[], it remains an interesting open problem 
hnding the deformations of the geometry which correspond to the addition 
of fundamental matter to the gauge theory. We have outlined a possible way 
to describe the required geometrical transition. We plan to return to this in 
a future work. 

Appendix 

In this appendix, using heuristic arguments, we propose a computational 
shortcut in the matrix model which yields the gauge theory effective super¬ 
potential. 

The proposal implies that the effective superpotential is written as a sum 
of terms, each arising from one particular held in the theory. We will be 
concerned only with the contribution of matter helds. As in the usual DV 
approach, the gauge held dynamics is taken into account by the addition 
of the Veneziano-Yankielowicz superpotential Wvy to the matrix model free 
energy. 

Our proposal is that the free energy of the matrix model associated to 
some gauge theory is given by the logarithm of its non-normalized partition 
function. Then, the gauge theory ehective superpotential is given by the 
(extended) DV relation: 

WeffecUve = WvY + -^ + y^X=l ' ( 45 ) 

An immediate consequence of this proposal is that a held not appearing 
in the superpotential contributes a factor of its “integration volume”: 

J DTi = Vol(Ti) . (46) 

As usual, we identify the’t Hooft coupling in the matrix model, g^M,[where 
M is dimension of the matrices) with the glueball superheld S. 
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In the case of massless adjoints, their volume of integration is the volume 
of SU{M) multiplied by that of a noncompact piece. When properly regu¬ 
larized the volume of the noncompact piece becomes an arbitrary scale. This 
contribution comes in the numerator, and can cancel the “volume of gauge 
group” which gives the VY term. If the gauge theory was N = 2 supersym¬ 
metric, there should not be any superpotentialQ for the glueball superheld S. 
This cancellation can be achieved if one identihes the cutoff in the volume of 
the adjoint held with the dynamically generated scale, A. 

Let us consider now the case of massless helds in the fundamental repre¬ 
sentation. Their integration volume is also inhnite, and one can also handle 
it by imposing a cutoff. However, it is not a priori clear what this cutoff 
should be identihed with in the gauge theory. It is nevertheless clear that 
one should use an object invariant under all symmetries of the (gauge) the¬ 
ory which also has the right dimension. A natural choice seems to be the 
determinant of the meson held: 

detX = |QQU,, . (47) 

To compute the integration volume for one quark we notice that the 
measure dQ dQ spans two fundamental representation; their product can be 
thought of as giving an adjoint representation and a singlet. The volume of 
the adjoint is hnite, but the volume of the singlet is inhnite. If we regularize 
this volume by imposing a cutoh at a “radius” X , the total volume will just 
be volume of the adjoint multiplied by ^ 

For each of the Nf quarks not appearing in the superpotential, one has to 
consider the corresponding integration volume. The cutoh for each volume 
will be interpreted in the gauge theory as the meson corresponding to that 
quark. Due to the various global symmetries of the theory these cutohs 
should combine to form an invariant, e.g. det X. It is not hard to see that, the 
net result of these considerations is a Weff{S,Xij) which yields the Aflleck- 
Dine-Seiberg superpotential upon integrating out the glueball superheld. 
Before integrating out the glueball superheld, the superpotential can be 

®The full analysis of the N = 2 dynamics is far beyond the scope of this proposal, see 
for steps in that directions. 

®This can be easily understood by analogy with the volume of a ball. The ball can be 
thought as a product of a sphere and i?+. The total volume is the product of the volume 
of the sphere and the radius to a power equal to the dimension of the space. 
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justified to be given by: 


W = NrS 



-NfS 



( 48 ) 


Indeed, the first term is the contribution of the gauge held; the second term 
comes from the compact part of the integration volume of the quarks; hnally, 
the last therm is given by the non-compact directions in the quark integration 
volumes, combined in a way compatible with the global symmetries of the 
theory. 

The mesons Xi appear now as cutoffs (regulators) in the integration vol¬ 
umes, and each gives a contribution to JF of the form logXj, as explained 
above. This justihes the X dependence of (^Sj) , as the superpotential is 


related to the hrst derivative of X with respect to S, 

When Nf = Nc, the hrst two terms cancel, and the only part left is 
Integrating out S gives a vanishing superpotential as well as the 


detX 


*51og^^. 
constraint 


=detX 


(49) 


which is the correct gauge theory result in the absence of baryons. 
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